o 
o 

(N 

-t— > 
o 

O 

o 

(N 

(N 
> 

in 

(N 

^: 
o. 
^• 

O; 

Oh! 
I 

Oh; 



X 



Comparative study of radio pulses from simulated hadron-, electron-, and 
neutrino-initiated showers in ice in the GeV-PeV range 

Shahid Hussaiifl and Douglas W. McKaj|t| 
Department of Physics & Astronomy University of Kansas, Lawrence, KS 66045. 

(Dated: February 2, 2008) 

High energy particle showers produce coherent Cherenkov radio emission in dense, radio- 
transparent media such as cold ice. Using PYTHIA and GEANT simulation tools, we make a 
comparative study among electromagnetic (EM) and hadronic showers initiated by single particles 
and neutrino showers initiated by multiple particles produced at the neutrino-nucleon event vertex. 
We include all the physics processes and do a complete 3-D simulation up to 100 TeV for all showers 
and to 1 PeV for electron and neutrino induced showers. We calculate the radio pulses for energies 
between 100 GeV and 1 PeV and find hadron showers, and consequently neutrino showers, are not 
as efficient below 1 PeV at producing radio pulses as the electromagnetic showers. The agreement 
improves as energy increases, however, and by a PeV and above the difference disappears. By 
looking at the 3-D structure of the showers in time, we show that the hadronic showers are not 
as compact as the EM showers and hence the radiation is not as coherent as EM shower emission 
at the same frequency. We show that the ratio of emitted pulse strength to shower tracklength is 
a function only of a single, coherence parameter, independent of species and energy of initiating 
particle. 

PACS numbers: 96.40.Pq,95.85.Bh,95.85.Ry,29.40.-n 



I. INTRODUCTION 



Ultrahigh energy (UHE) neutrinos, neutrinos with en- 
ergies on the order of a PeV and above, are predicted 
at some flux level by most models of the highest energy 
astrophysical processes observed in the Universe. More- 
over, observations of cosmic rays above energies vifhere 
gamma-nucleon thresholds open up imply the existence 
of cosmogenic, or Greisen-Zatsepin-Kuzmin (GZK), UHE 
neutrinos[l|. Ultrahigh energy neutrino detection will 
probe interactions at energies well beyond the reach 
of man-made accelerators, while bringing us directional 
and internal information about the fantastically energetic 
sources. 

A standard method of UHE neutrino detection relies 
on showers initiated by the electrons and hadrons pro- 
duced in neutrino-nucleon interactions. Unlike muons 
and taus, electrons and hadrons dump their energy very 
quickly in matter through electromagnetic and strong in- 
teractions. These showers produce radio and optical sig- 
nals 2] , can be detected in detectors like RICE 3, 4] and 
AMANDA 0, and by those near deployment like ICE- 
CUBE 111, under development like ANITA 0, or in out- 
line stage like SALSA % or X-RICE '^. The high degree 
of transparency of cold ice and salt to radio transmission 
make Antarctica or large salt domes attractive for radio 
detection schemes like RICE, ANITA or SALSA. This 
paper extends earlier studies on the signal production 
and propagation in cold ice. It is readily adaptable to 
salt, as well. 



Though considerable work has been done on shower 
simulations of single particles in ice in different approx- 
imations at high energies jlj, uMt uM ' most of the effort 
has focussed on EM showers and the pulses they produce. 
The initial flux limit from the RICE group, for example 
PI, reported an upper limit on i/g flux alone. This was 
conservative, in the sense that only the electron shower 
was used in calculating the pulse from a charged current 
(CC) interaction, without the hadronic component and 
with no neutral current (NC) interactions of Ve included. 
No ly^ or Vr interactions were treated in the main analy- 
sis, though estimates of hadronic contributions were dis- 
cussed and presented in an auxiliary table. More recently 
the preliminary results of the analysis on an expanded 
subset of data |4| were reported, in which the hadronic 
contributions to the pulses from Vg CC were estimated 
and included in a preliminary, updated limit. Detailed 
hadronic shower and pulse studies and determination of 
the degree to which hadronic and electromagnetic com- 
ponents are coherent were not available in that study. 
Toward filling this gap we report here: (i) Full 3-D sim- 
ulations of EM showers to a PeV and hadronic showers 
to 90 TeV for single hadrons. (ii) Direct, full simulation 
of the showers initiated by charged current interactions 
of all neutrino flavors, up to a PeV. (iii) Time struc- 
ture analysis of the EM, hadronic, and neutrino showers 
and discussion of the sources of difference between their 
pulse strengths, (iv) A universal relationship between 
track length and emitted pulse strength, independent of 
shower type or energy. 
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II. OVERVIEW OF PYTHIA/GEANT 

SIMULATIONS AND CALCULATION OF 

SIGNAL STRENGTH 

For present purposes, we take our target and detection 
volume to be a kilometers-cubed volume of cold ice, with 
an instrumented volume that is a small fraction of this. 
Operationally, the calculations break into two stages: (1) 
neutrino interactions with hadronic targets to produce 
a lepton or neutrino plus multi-hadron final state (with 
PYTHIA _13]) (2) an electromagnetic (EM) and hadronic 
shower (from Vf^ charged current (CC) events) or only 
hadronic shower (from all neutral current (NC) events 
and v^, and z/^ CC events jl^) (with GEANT [3) and 
the real time calculation of the EM pulse generated by 
the shower particles. The output field can be described 
as a pulse of field a few nanoseconds long moving out 
from the interaction point on a Cherenkov cone of half 
angle 56*^ and Gaussian width of several degrees [lO|. 
The width of the cone is frequency dependent and, above 
the Landau-Pomeranchuck-Migdal (LPM) effect thresh- 
old, energy dependent IL]. The LPM threshold is above 
the energies we cover here. 

To explore hadronic shower-generated pulses in detail, 
we expand considerably on earlier GEANT- based shower 
and pulse studies [l3|. We use Geant4.5.2 to push EM 
shower analysis to 1 PeV and also include the same anal- 
ysis of TT''" , 7r~ , and proton hadron-induced showers near 
100 TeV. Beyond 100 TeV, the energy limitations on 
some of the input hadron cross sections prevent full de- 
velopment of the shower with GEANT4.5.2 [ll |. 

To simulate neutrino induced showers and the radio 
frequency pulses they produce, we use PYTHIA6.205 to 
simulate neutrino-nucleon CC interaction. The PYTHIA 
output data files are fed to GEANT and relevant shower 
and field parameters calculated as the shower simulation 
proceeds. The Fourier transform of the electric field in 
the Fraunhoffer limit is calculated for each track segment 
of every particle in the shower and the contributions of 
the track segments to the total field are summed vec- 
torially, as described in [l^. We follow the standard 
practice and report the magnitude of the field at a given 
angle and frequency multiplied by the distance from the 
start position of the shower to the position of the ob- 
server. This value is a function of frequency, angular po- 
sition of the observer and energy of the initiating particle. 
We generally refer to this as "the pulse" . The phase of 
the field produced by each track segment depends upon 
the track start and stop times and initial track position, 
which feeds local information into the total pulse. 

We report here the essential features of showers and 
the pulses they generate. We want to find out whether 
the features of hadron induced showers and pulses at a 
given energy are similar to EM features at the same en- 
ergy, whether the comparison is energy dependent, and 
whether the multi-hadron final states from z/e CC inter- 
actions produce pulses that are essentially coherent with 
the pulse the electron produces. We consider electron 




Log(Y) = A + B • Log(X), 
A (intercept) = 3.129 +/- 0.002 
B (slope) = 0.998+/- 0.001 
R (correlation coefficient)- 1 



E(TeV) 



FIG. 1; Projected weighted-track length difference r~ — rf 
(see text for more detail) vs energy for different particles. 
Equation of the linear fit for e~ is also shown. For the neu- 
trino case, elasticity (1 — y) is 0.65+/-0.01. 



and i>e CC induced events first and then turn to hadron 
and ly^j and Vr CC-induced events. The simulation re- 
sults of CC-induced events of f^ and i^t can also be used 
to estimate the results for NC-induced events as both are 
initiated hadronically (tau and muon decay lengths are 
large as compared to the typical shower depths at high 
energies. Hence, for a given value of inelasticity, we ex- 
pect Vf^ and i^r CC-induced showers to be scalable to the 
NC-induced events with the same inelasticity value). 

The simulations were done in 3-D, including all the 
hadronic and EM physics processes involved. For each 
particle except neutrinos, we have simulated 100 show- 
ers at 100 GeV, 50 showers at 1 TeV, 10 showers at 10 
TeV, single showers at 90/100 TeV and 1 PeV. For the 
neutrinos we have simulated 100 showers at 1 TeV, 30 
showers at 10 TeV, and 5 showers at 100 TeV. The mean 
deviations of the non-EM numbers given in Tableware in 
the range 1-8%; although we cannot directly determine 
the error in the data at 90 TeV and above, as we have 
only single showers, we expect it will certainly be below 
10% because a single shower will have many times the 
number of particles in a 10 TeV shower, so the relative 
fluctuations will be smaller in a shower around 100 TeV 
and above. 



Electron and Electron-Neutrino Charged 
Current Induced Pulses 



As in earlier work all field calculations are done in the 
Fraunhoffer limit i? » -^, where R is detector-shower 
distance, A is the wavelength, and L is the shower pene- 
tration distance. The simulation is run in cold ice with 
p = Omg/cm? , and with A = 14.3, Z = 7.2, and radia- 
tion length=38.8cm, as calculated by GEANT. 

Injecting electrons at 0.1, 1, 10, 100, and 1000 TeV, we 



TABLE I: Different characteristics of the showers as a func- 
tion of energy; 1st column is energy in TeV; 2nd column has 
the particle name (for neutrinos the elasticity value is 0.65 in 
every case) ; Ar^ — r~ —rf in the 3rd column gives the differ- 
ence of the track lengths in meters, for positive and negative 
charges, projected along z direction (i.e. direction of prop- 
agation of the primary particle); 4th column (_R|_E|1) gives 
R\E\ (see text for definition of -R|-E|) values in iiV/MHz at 
the Cherenkov angle, at 1 GHz (see Section II for the es- 
timated random errors of these numbers); 5th column gives 
w = -^f- (its meaning is discussed in section III); 6th column 
gives b{w) * lO"" (see Eq. 2); 7th column gives R\E\2 which is 
the pulse value calculated from Eq. 2. 



E 


name 


Ar^ 


R\E\l 


w 


b{w) * le5 


R\E\2 




e~ 


136 


0.014 


0.90 


9.9 


0.013 


0.1 


TT~ 


70 


0.0034 


0.62 


4.7 


0.0033 


t:+ 


68 


0.0034 


0.62 


4.7 


0.0032 




P 


53 


0.0026 


0.56 


4.8 


0.0026 




e~ 


1336 


0.13 


0.90 


9.9 


0.13 




1T~ 


898 


0.058 


0.75 


6.4 


0.058 




t:+ 


897 


0.057 


0.74 


6.3 


0.056 


1 


P 


812 


0.046 


0.71 


5.7 


0.046 




l^e 


1156 


0.087 


0.81 


7.8 


0.091 




'^M 


269 


0.012 


0.60 


4.6 


0.012 




Vr 


271 


0.012 


0.60 


4.6 


0.012 




e~ 


13319 


1.3 


0.90 


9.9 


1.3 




TT~ 


9963 


0.75 


0.79 


7.4 


0.73 




7r+ 


9728 


0.70 


0.78 


7.1 


0.69 


10 


p 


9210 


0.64 


0.77 


6.9 


0.64 




Ve 


12095 


1.1 


0.86 


9.0 


1.1 




^M 


3296 


0.22 


0.77 


6.9 


0.23 




Vt 


3395 


0.25 


0.78 


7.1 


0.24 


10^ 


e^ 


1.33S5 


13 


0.90 


9.9 


13 


90 


■K~ 


1.04B5 


9.2 


0.85 


8.8 


9.1 


10" 


t:+ 


1.19S5 


11 


0.86 


9.0 


11 


90 


P 


90007 


6.9 


0.80 


7.5 


6.8 


10^ 


Ve 


1.26S5 


12 


0.87 


9.3 


12 


10^ 


^M 


39982 


3.5 


0.84 


8.5 


3.4 


10^ 


Vt 


36800 


3.0 


0.82 


8.0 


3.0 




e~ 


1.33S6 


132 


0.90 


9.9 


132 


10^ 


Ve 


1.32S6 


131 


0.89 


9.7 


128 


Uf, 


4.08B5 


37 


0.86 


9.0 


37 




Vr 


3.67B5 


30 


0.82 


8.1 


30 



calculate Arz = r~ —r'^ , the total weighted track lengths 
projected along shower axis, shower charged particle pro- 
files, frequency dependence of the field {\E\) times the 
distance R to the detector (the field strength decreases 
as 1/R in the Fraunhofi^er limit), and angular dependence 
of the field at 1 GHz. The results for key parameters are 
shown in Table Q] . The electron is the first entry at 
each energy. As is apparent from the numbers, the pro- 
jected track length difference (Ar^) between the negative 
and positive particles and the pulse height at a fixed fre- 
quency and at the Cherenkov angle grow linearly with 
energy to high accuracy. This is shown in Figs. ^ and |2 
where the open boxes are the electron data and the fits 
are shown on the graph. These results agree with all ear- 
lier work lid Il2| . The angular distribution is shown in 




fit for e-: 

Log(Y) = A + B • Log(X), 

A (inlercept) = -0.864 +/- 0.006 

B (slope) = 0.993 +/- 0.003 

R (coiTelalion coel'ridenl)= 0.9999S 



E(TeV) 



FIG. 2: R\E\ aX 6c and 1 GHz vs energy for different parti- 
cles. Equation of the linear fit for e~ is also shown. For the 
neutrino case, the elasticity is the same as in figure 



Fig. O for E = IPeV in the region of the Cherenkov an- 
gle. This plot clearly shows the rapid suppression of sig- 
nal off the Cherenkov cone. Injecting 1, 10, 100, and 1000 
TeV neutrino-nucleon interaction output from PYTHIA, 
we calculate the corresponding quantities with elasticity 
I — y = 0.65 ± 0.01 at each energy. We choose the same 
y at each energy so that we can determine whether the 
hadronic component of the entire shower becomes more 
efficient at generating pulse strength with increasing en- 
ergy. This requires keeping the hadronic fraction of the 
total energy constant as we change energy. The choice 
y = 0.35 turns out to be big enough so that hadronic 
fluctuations from shower-to-shower are reasonably small 
and at the same time the electron does not completely 
dominate the signal in the Ve case. The trends in the 
behaviour with neutrino energy that we describe below 
are true regardless of the y- value chosen. 

In Table m, selected examples of these calculations are 
given at each energy. As is well known 15], the elas- 
ticity peaks at 1 and the average grows slowly with en- 
ergy, starting at ???? at 0.1 TeV and reaching about 
0.74 at Ei^ — IPeV. In incident neutrinos of given fla- 
vor and energy, we assembled a sample within our chosen 
range of y by generating PYTHIA events and selecting 
those in our chosen y range until our sample size was 
reached. As mentioned above, our aim is to study the 
pulse-production efficiency of the hadronic component of 
neutrino showers as a function of energy. In other words, 
is the pulse growth slowers than, the same as, or faster 
than the growth in hadronic energy? We need the frac- 
tion (y), of hadronic energy in the shower to be held 
constant as we raise the energy to answer this question. 
Since average y decreases as the energy increases, the av- 
erage y at a given energy is not the best choice for our 
purpose here. For application to simulation of signals 
for an actual experiment, one would not use fixed y of 
course. At a given y value, the fiuctuations in shower and 




FIG. 3: R\E\vse,Sitl GHz, and E=l PeV for single shower 
simulations. For neutrinos, the elasticity is the same as in 
figure 



pulse parameters decrease with increasing energy and at 
the same time the Ve results look more and more like 
electron results; the Ve event at 1 PeV gives projected, 
weighted track length and pulse values that are essen- 
tially the same as those of the electron event. This is 
clearly shown in Fig. |3as well. 

The trend of our simulation data shows that, as the 
primary energy increases, the hadronic component of the 
events for any y value contributes in such a way that the 
Vf, events look more and more like electron events as the 
primary energy increases. We will see this in a more di- 
rect way in the results from single injected hadron, which 
we describe in the next subsection. 



B. Hadron, and u^ and i/t CC Induced Pulses 

The calculation of the charged pion and proton induced 
events proceed similarly to those of the electron/positron 
induced events. There are many more processes in- 
cluded in the hadronic cases, of course; the maximum 
energy of pion and proton showers is around 100 TeV 
in GEANT4.5.2, because some cross sections, such as 
the kaon- total cross sections, are extrapolated only to 
about 10 TeV. Nonetheless, we can gain valuable in- 
sight from the comparison between electron and hadron- 
induced events. Since the low energy tail of a shower is 
essentially EM, and the excess charge develops in this low 
energy region, the hadron shower should behave similarly 
to the EM showers to the extent that gammas from 7r° 
decay dominate the end of the shower. This should be 
a better picture as energy increases and fluctuations in 
particle species wash out. 

The essential features of the comparisons between the 
electron and hadron-induced showers are summarized in 
Figs. Hand [3 The projected weighted track length and 
pulse height at Cherenkov angle at fixed frequency grow 



linearly with energy. It is clear from the graphs (Figs. Q] 
and |2Jl and from a look at the numbers themselves (Ta- 
ble P) that the efficiency of converting hadron input en- 
ergy increases with energy, as we emphasize further be- 
low. It is important to note that a given weighted pro- 
jected track length in a hadronic shower is not as efficient 
in producing the radio signal as is the same length of EM 
shower. However, this difference between the hadronic 
and the EM showers shrinks at higher energies. As we 
will discuss in the next section, this is because of the 
difference between hadronic and EM showers in the dis- 
tribution in space of the particles contributing to the 
weighted, projected track length. 

As we see in Fig. |21 the pulse outputs from all the 
showers are nearly equal as one approaches 1 PeV shower 
energy. We expect that the hadron curve will flatten out 
and follow the electron at higher energies, though we can- 
not verify this with the current GEANT capabilities. The 
tendency for hadron events to act more like EM events is 
consistent with our observation earlier that the Ve events 
act more and more like electron events as we raise the 
injection energy; the hadron component adds coherently 
with the EM component. 

The Vf^ and i>r CC results are the final entries at each 
energy. Here again it is clear that the hadronic com- 
ponent of the final state is becoming more efficient in 
converting energy to EM pulse strength as the injection 
energy increases. For example, the average v^^ or i/,- event 
at 10 TeV with our illustrative value \ — y = 0.65 pro- 
duces a significantly larger fraction of the electron pulse 
than it does at 1 TeV. The corresponding comparison at 
1 PeV for Vi^ and Vr events produce noticeably different 
shower and pulse results at the same energy and elastic- 
ity. We looked at the detailed structure of the PYTHIA 
output for these and lower energy cases and found that 
the larger track length and larger pulse result correspond 
to the event with larger 7r° content in the PYTHIA out- 
put (GEANT4 input) (and, ultimately, higher shower 
EM-content; this explains why v^^ produces a stronger 
pulse than v^ with the same E and y in one example and 
weaker in another. Even at 1 PeV, there are significant 
fluctuations in the species and subenergy compositions of 
the u + p ^r l-\- X final states and these fluctuations pro- 
duce corresponding shower and pulse fluctuations. The 
trends in Figs. ^ and |21 clearly indicate that above 1 PeV 
these fluctuations damp out and hadronic components of 
the showers act like electromagnetic showers of the same 
energy. The LPM effect is expected to set in at higher en- 
ergies and elongate the EM showers but not the hadronic 
showers [ll|. 

We have also calculated the frequency spectrum of the 
pulses at Cherenkov angle from each type of shower at 
each energy. In summary, we find the spectra from differ- 
ent types of shower look more and more similar to each 
other as the shower energies increase. By 10 TeV, they 
all rise sharply from zero at zero frequency, peak at sev- 
eral GHz, fall and then level out in the 10 to 20 GHz 
region and are essentially flat above that. The frequency 



behavior up to 1 TeV is discussed for electron showers in 
S. Razzaque et al. [I2|, where, as in our present study, 
the sharp turn over of the spectrum does not set in until 
1 TeV. We find that the hadron and hadron dominated 
v^ and Vr do not show this behavior until and order of 
magnitude higher energy. 



comparison. 

Finally we examined the degree to which the value of 
r J — r+ is determined by the particles which are highly 
relativistic and are hugging a point moving with the 
speed of light along Z-axis. For energies of 0.1 TeV, 1 
TeV, 10 TeV and 90 TeV [l3, we compute the value of 
r~ — r+ for just those tracks with points within a distance 



III. TIME STRUCTURE ANALYSIS OF 
HADRONIC AND EM SHOWERS 

As far as we are aware, direct time development of EM, 
hadronic and neutrino-induced showers in ice has not 
been reported. Here we describe several analyses to clar- 
ify the reasons behind the generally improved efficiency 
in radio pulse production by hadronic showers of increas- 
ing energy. We focus on several direct-time evolution 
representatives of the shower content for this purpose. 



= v/^2 + Y^ + {Z - cty < 0.15TO, 



(1) 



The weighted, projected track length Ar^ —r 



de- 



scribed earlier, is known to be a good indicator of pulse 
height produced by the EM showers ^3| • As Table U 
shows, this is not the whole story in understanding the 
difference between the pulses produced by electron show- 
ers vs. hadron showers at lower energies. For example. 



at 1 TeV, the ratios of r^ 



for e : TT , e : 7r+, and 



e : p are 1.49, 1.49, and 1.65, respectively, while the 
corresponding R\E\ ratios are 2.24, 2.48, and 2.83. The 
corresponding comparisons for Ve, Vfi and Vt are 1.16, 
4.95 and 4.95 compared to 1.53, 11.0 and 11.0. The ra- 
tio of pulse to track length is larger in the EM than in 
the hadronic showers and, to the extent they are domi- 
nated by hadrons, the neutrino showers. The EM showers 
are evidently producing more coherence in the Cherenkov 
emission. Intuitively one expects the degree to which the 
emission is coherent depends upon the degree to which 
the charged particles, and therefore the charge imbal- 
ance, are concentrated in a small volume, acting like a 
single, point charge. 

To get a visualization of these features as the shower 
develops, in Figs. 4 and 5 we plot the X, Y, and Z compo- 
nents of the charged particle position versus time at each 
step in every track's development in a 100 GeV electron 
and proton showers. The lightest shaded points are X, 
next darker Y, and the darkest is the Z component of all 
charged particle positions as a function of time. Beyond 
50 ns, there are no tracks in electron showers and only 
a tiny percentage in hadron showers. The speed-of-light 
line in cold ice is designated by a heavy dark line. A 
comparison between the plots vividly shows that the X 
and Y excursions of electron shower particles are smaller 
than those of the proton shower of the same energy. The 
concentration of the Z components near the speed of light 
in vacuum is correspondingly higher than that in the pro- 
ton shower. These plots indicate that the coherence of 
the emission from the electron shower should be greater 
than the one from a proton shower, which correlates well 
with the fact that the emission is weaker from the proton 
shower at 1 GHz than is expected just from the rj — r+ 



as shown in Fig. 6 (only two dimensions are shown). The 
choice a = 0.15m is guided by the wavelength at 1 GHz. 
This choice, though not unique by any means, is satisfac- 
tory to make our point about the role of compactness and 
coherence |23- The idea is simply that charge dispersed 
farther from the shower core will tend to degrade coher- 
ence. We have calculated the projected weighted-track 
lengths Ar^ = r~ —r^ of the particles within the shaded 
region of the sphere with radius a=0.15m (see figure 6) 
and the projected weighted-track lengths Ar^ = t~ — r'^ 
due to all the charged particles in the shower. 

In Table m we show the ratio w = -jr^ for the showers 
initiated by electrons, pions, protons and the three flavors 
of neutrinos at various energies. As we see here, this ra- 
tio has a constant value of 0.90 for the electron initiated 
showers between energies 100 GeV-1 PeV. This means 
90% of the projected weighted-track length comes from a 
compact region defined by Eq. ^ On the other hand this 
ratio is smaller for a proton, a pion, or a neutrino shower 
and it improves as we go to higher energies. For EM 
showers, the value of Arx has long been recognized as an 
extremely reliable measure of the pulse height produced 
at Cherenkov angle for frequencies at and below 1 GHz. 
We emphasize here the added insight gained by treating 
separately the track length of the core region defined in 
Eq. 1, and the track length outside this region. If we 
look at Tableland the Figs. 2 and 3, we see as the ratio 
w for hadronic showers approaches that of the electron 
shower (TableQ]), the weighted, projected track lengths 
and pulse heights of the hadronic showers also approach 
those of the electron showers. The well established lin- 
ear relationship between electron shower projected track 
length and energy and between pulse height and energy is 
expressed in terms of our "compactness" parameter as a 
constant value w — 0.90 for all electron showers between 
0.1 TeV and 1.0 PeV. 

As is clear from the numerical results summarized in 
Tableland our discussion above, the pulse strength at 
the Cherenkov angle at a given frequency is related to 
the weighted, projected track length, Ar^, and the ratio 
w. The data suggest that the ratio of R\E\ to Ar^ is 
a function of w that is independent of the particle type 
and energy. For example, there are three entries with 
w = 0.86: i/e at 10 TeV, 7r+ at 100 GeV, and v^ at 1 
PeV; these have values i?|£:|/Ar^ = 9.1 10-^ 9.2 10"^ 
and 9.1 10~^, respectively. This and similar results in 



Table U indicate there is a relationship 



R\E\/Ar,^b{w), 



(2) 



where w is defined above and one determines the func- 
tion b{w) from the data. We find that a third degree 
polynomial in w fitted to the data plotted as R\E\/Arz 
vs. w gives a remarkably good description of the data in 
Table Q] |2j| . Given the value of w and Ar^ , one repro- 
duces the value of i?|-E| within a few percent for all but 
a couple of the 29 entries, and those two are good to bet- 
ter than about 6%. As remarked in footnote 21, above 
1 TeV, a linear fit works quite well. This universal fea- 
ture supports our picture that the key factor that relates 
the pulse at a given frequency to the weighted, projected 
track length for showers induced by any particle in the 
energy range we consider is the degree of coherence of 
the emission, as determined by the compactness of the 
relativistic core of the shower as it moves along through 
the medium. 



IV. SUMMARY AND CONCLUSIONS 

Ultrahigh energy neutrinos will provide a uniquely 
clear window to explore fundamental physics and astro- 
physics processes at the highest known energies. These 
neutrinos produce huge particle showers when they in- 
teract with dense matter. The showers then produce co- 
herent, Cherenkov radio signals, the strength of which 
depends on the net charge in the shower, on the track 
length projected along the shower axis and on the de- 
gree of compactness of the shower. All of these factors 
increase with neutrino energy. Combined with a radio- 
transparent medium such as cold ice or pure salt, the 
radio detection technique becomes increasingly powerful 
as energy rises in the UHE regime. 

The particle content early in a shower depends on 
the flavor of the neutrino producing it. Electron- 
neutrino showers immediately have a heavily electromag- 
netic character, which dominates the hadronic compo- 
nent on average. On the other hand, the muon-neutrino 
and taon-neutrino showers are essentially hadronic in 
character. Below a PeV our simulation results show that 
hadronic showers are not as efficient in producing radio 
signals as EM showers. Though the pulse strengths pro- 
duced by hadron showers increase with energy, the net 
charge, the weighted, projected track length, and the cor- 
responding pulse strength are all smaller than those of 
the EM showers at the same energy. In addition, we 
found that the pulse strength per weighted, projected 
track length is smaller for the hadronic showers than for 
the electromagnetic ones. However, the key result of our 
study is that the differences in the efficiencies of pulse 
production disappear as one goes to higher and higher 



energy. The same is true of the comparisons between the 
neutrino showers and the electron showers. Extrapolat- 
ing hadronic shower results to 1 PeV from our highest 
simulation results at 100 TeV, we find that the hadron 
and electron showers at a PeV give essentially the same 
pulse strengths. Moreover, at 1 PeV, the simulated v^ 
shower with elasticity of 0.65 is indistinguishable from 
the electron shower at the same energy. The i^^ and i>r 
showers continue to show convergence toward electron 
shower behaviour when the elasticity factor in the events 
is accounted for. 

In a new type of study employing time- and space- 
structure of the showers, we showed that the low energy 
electron showers are significantly more compact than the 
corresponding hadronic showers and neutrino showers at 
low energies. The coherence is directly related to the 
compactness and therefore to the efficiency of converting 
charge imbalance and track length into EM pulses. We 
show that a simple, universal formula nicely describes 
this aspect of the simulation data. Again, as one goes 
to ultrahigh energies, hadronic and EM showers become 
equally compact and equally coherent and, consequently, 
equally efficient in producing radio signals . All the sig- 
nals rises linearly with energy at ultrahigh energies and 
with essentially the same normalization. We believe our 
study develops a detailed account not previously avail- 
able of how his simple picture emerges in the 100 GeV to 
1 PeV energy regime. 

We conclude that at and above a PeV, until 
the Landau-Pomeranchuk-Migdal (LPM) effect becomes 
important [Ig, the hadronic and electromagnetic showers 
produce essentially the same pulses at a given energy. 
In addition, the hadronic and electromagnetic compo- 
nents of the Ve pulses are effectively coherent. Because 
the LPM effect has bearing on the EM particles that are 
above LPM threshold, and the EM particles that evolve 
from UHE hadronic showers up to EeV energies are below 
this threshold [ij , the hadronic showers are predicted to 
become more effective than the EM showers at producing 
EM pulses at some point above the LPM threshold. Fur- 
ther study of this important problem is certainly needed 
and will be the subject of future research. 
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Figure Captions for Figures 3-6 
(separate .jpg files) 

FIG. 4: 3-D structure in time of a lOOGeV electron 
shower. The light gray color represents X, dark gray Y, 
and black Z coordinates. For the Z case, particles below 
the thick solid line z=c*time/1.8 are moving slower than 
the speed of light in cold ice. 

FIG. 5: Same as figure 4 for a lOOGeV proton shower. 

FIG. 6: The circle shown here has a radius a=0.15m 
and is centered at a point that moves with the speed of 
light in vacuum. The shaded region is defined by Eq. 1 
in the text. 
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